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Abstract: Previous research on gene expression analysis and association tests have suggested that RELN is a risk 
gene for schizophrenia in world populations. Based on the reported down-regulation of RELN in schizophrenia patients 
compared with normal subjects, we speculated that variants in the RELN promoter region may confer risk for 
schizophrenia. In this study, we investigated the associations of three SNPs in the promoter region of RELN with 
schizophrenia in a case-control sample from southwestern China (940 cases and 1369 controls). The results suggested that 
none of the SNPs showed significant associations in our sample, indicating the risk variants for schizophrenia in RELN 
may not be located in the promoter region. We also performed meta-analysis by combining our data with previously 
reported data on the Chinese population with a total sample size of 2843 individuals, and the result remained non- 
significant. Collectively, our results suggested variants in the RELN promoter may not harbor risk SNPs associated with 


schizophrenia in the Chinese population. 
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Schizophrenia is a severe psychiatric disorder with et al, 2007; Williams et al, 2004; O'Donovan et al, 2008). 


high heritability. Many susceptibility genes have been However, many cannot be successfully replicated among 
proposed by linkage analyses, candidate gene studies, different populations. A recent GWA study on 
and genome wide association (GWA) analyses, including Ashkenazi Jews identified that a common variant 


DISC1, NRGI, TCF4, NRGN, and ZNF804A (Thomson (rs7341475) of RELN was significantly associated with 
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schizophrenia in females (Shifman et al, 2008). This 
association was further supported by replications in a UK 
study and another Ashkenazi Jew study (Shifman et al, 
2008; Liu et al, 2010), but was unable to be replicated in 
USA and Chinese populations (Shifman et al, 2008; 
Need et al, 2009), 
heterogeneity in RELN for schizophrenia between 


indicating potential genetic 
different world populations. 

We previously identified several RELN SNPs and 
haplotypes associated with schizophrenia in the Chinese 
population (Li et al, 2011), which is consistent with 
reported associations in European populations and 
suggests RELN is likely a common risk gene for 
schizophrenia in populations worldwide, though the risk 
variants differ between reported associations (Kahler et 
al, 2008; Wedenoja et al, 2008, 2009). 

Down-regulation of RELN mRNA and protein in the 
brain of schizophrenia patients has been reported 
previously, but the underlying mechanism remains 
unknown (Impagnatiello et al, 1998). Hypermethylation 
of RELN promoter in schizophrenic patients is a likely 
cause; however, negative findings have also been 
reported (Grayson et al, 2005; Tochigi et al, 2008). One 
reasonable explanation is that mutations in the RELN 
promoter may affect gene expression and function and 
influence neurodevelopment, leading to schizophrenia 
susceptibility. Interestingly, short tandem repeats in the 
promoter region influence RELN promoter activity, but 
the association tests in schizophrenia show negative 
(Akahane et al, 2002). Additionally, an 
investigated RELN promoter SNP has shown marginal 


results 


significance with schizophrenia in a small case-control 
Chinese sample (P=0.08) (Chen et al, 2002), suggesting 
the RELN promoter may harbor risk variants for 
schizophrenia. In this study, we examined the association 
between the RELN promoter SNPs and schizophrenia in 
a Chinese case-control sample. 


1 Methods and Materials 


1.1 Samples 

We recruited 940 unrelated schizophrenia patients 
(481 females and 459 males, mean age=37.4 a, SD=9.1) 
and 1 369 unrelated normal controls (732 females and 
637 males, mean age=36.7 a, SD=6.8) from southwestern 
China. The patients were all from Yunnan Mental Health 
Hospital and The Second People's Hospital of Yuxi City 
and were diagnosed with having schizophrenia according 
to DSM-IV and ICD-10 criterion. Patients who had 


history of alcoholism, substance induced psychotic 
disorders, epilepsy, neurological diseases, or other 
symptomatic psychoses were excluded from this study. 
Control subjects were recruited from the local general 
populations. All individuals were asked to provide 
detailed about 
psychiatric histories, and those who had a history of 


information medical and family 
mental disorders, drug abuse, or alcohol dependence 
were excluded. All patient and control subjects were of 
Han Chinese origin from the Yunnan province of 
southwestern China. All individuals were provided with 
written informed consents for participation, and the 
research protocol was approved by the internal review 
board of the Kunming Institute of Zoology, Chinese 
Academy of Sciences. 

1.2 SNP selection 

As no SNPs were shown in the -2 kb region 
upstream from the transcription start site (TSS) of RELN 
in Han Chinese from Beijing (CHB) data obtained from 
HapMap database, we used bi-directional sequencing to 
search for potential SNPs in the promoter region of 
RELN on 100 schizophrenic and 100 healthy randomly 
selected subjects. Primers used for PCR amplification of 
the RELN promoter were 5'- AGCCCAGAAGCAATGAA 
TAAC-3' (forward) and 5'-TCCCAACTTGTGACTCCATT 
C-3'(reverse). The PCR program started with an initial 
incubation at 95 °C for 5 min, followed by 40 cycles of 
95 °C for 30 s, 60 °C for 40 sec and 72 °C for 1 min, and 
then held at 72 °C for 10 min. Two SNPs were identified 
(g.-847G>A and g.-888G>C), with g.-888G>C having 
been investigated previously (Chen et al, 2002). We 
further carried out case-control association analysis on 
these SNPs (plus another RELN promoter SNP 
1s6951875, shown in HapMap data) in our samples. 

1.3 SNP genotyping 

Venous blood was collected from all participants, 
and genomic DNA was extracted from the blood sample 
using the phenol-chloroform method. The DNA samples 
of the cases and controls were randomly distributed in 
the case-control DNA plates. 

All selected SNPs were genotyped by SNaPShot as 
described in our previous study (Luo et al. 2008). In brief, 
genomic fragments which contained selected SNPs were 
amplified by PCR with a total volume of 25 uL 
(including 10 ng of genomic DNA) in 96-well plates. 
Amplified fragments were purified and specific 
genotyping primers were used to amplify the target site. 
After one base extension, the reaction was terminated 
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and the products were loaded on an ABI 3130 automatic 
sequencer (Applied Biosystems). Primer sequences for 
SnaPShot analysis were 5'-TTTTTATGAGGTATTCTGA 
CACTGGATGAAGAATAATTAT-3'(ts695 1875), 5'-TTTT 
TITTTTTTTTTGCGGGAGGGACAGGGGGCCTGGGT 
GGGAAGGGAGC-3' (g.-847G>A) and 5'-TTTTTTTTTT 
TITTTTTTTTATGAGGCTCTGTCGCTGCCGCGAGGG 
GCCGGGCGG-3'(g.-888G>C). The SNP genotype 
callings were automatically performed using ABI 
GeneMapper 4.0 and verified manually. To ensure 
accuracy of genotyping, we used bi-directional 
sequencing on 100 randomly selected individuals. No 
genotyping errors were found. 
1.4 Statistical analysis 

The Haploview program was applied to test the 
genotypic distribution of SNPs for Hardy-Weinberg 
equilibrium (HWE) between paired SNPs, and to define 
haplotype blocks (Barrett et al, 2005). Allelic and 
genotypic associations were accessed with PLINK 
(Purcell et al, 2007). To detect significant differences in 
association in female or male samples separately, we 
conducted statistical analysis with sex as a covariate 
using PLINK (Purcell et al, 2007). The 95% confidence 
intervals (CI) of odds ratio were calculated with the 
(http://faculty. vassar.edu/lowry/odds2x2. 
html). Haplotype frequency estimation and association 


online tool 
tests were performed using PLINK, and only those 
haplotypes with a frequency of >0.01 in cases and 
control subjects were considered (Purcell et al, 2007). 
Power analysis was performed using G*power program 
(Erdfelder et al, 1996). For meta-analysis, we used the 
Mantel-Haenszel method with a fixed-effects model. 
Analysis was conducted by RevMan manager (The 
Cochrane Collaboration, 2002). 


Vol. 32 


2 Results 


Due to genotyping failure of some samples, analyses 
were based on 2 230 samples for SNP g.-847G>A (861 
cases and | 369 controls), 2 307 samples for SNP g.- 
888G>C (940 cases and 1 367 controls), and 2 296 
samples for SNP r1s6951875 (940 cases and 1 356 
controls). The overall genotype calling rate was 98.6%. 

Genotype distributions of the three SNPs in both 
cases and controls were in HWE (P>0.05). The LD map 
of the tested SNPs in cases and control samples are 
shown in Fig. 1. 
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Fig. 1 The LD map of the tested SNPs in case and 
control samples 
Linkage disequilibrium (LD) between the paired SNPs. The haplotype 
blocks were estimated by Haploview using the standard D’ confidence 


interval (CI) algorithm. 


None of the tested SNPs were significantly 
associated with schizophrenia (Tab.1 and Tab.2). We 
also investigated the association of these SNPs with 
schizophrenia in females and males separately, and no 


Tab. 1 Allele frequencies and single SNP association (P-value) analysis 











MAF 
Sample SNP Minor allele P-value OR(95%CI) 
Case Control 
g.-847G>A A 0.2189 0.2213 0.851 0.99(0.85-1.14) 
All samples g.-888G>C C 0.0803 0.0877 0.371 0.91 (0.73-1.12) 
rs6951875 T 0.0952 0.0867 0.319 1.11(0.90-1.36) 
g.-847G>A A 0.2184 0.2316 0.463 0.93(0.76-1.13) 
Female samples g.-888G>C C 0.0769 0.0840 0.532 0.91(0.67-1.23) 
rs6951875 T 0.0988 0.0815 0.144 1.24(0.93-1.64) 
g.-847G>A A 0.2195 0.2096 0.585 1.06(0.86-1.31) 
Mals g.-888G>C C 0.0839 0.0921 0.503 0.90(0.67-1.22) 
samples 

rs6951875 T 0.0915 0.0926 0.933 0.99(0.74-1.33) 


Abbreviations: MAF-minor allele frequency; OR-odds ratio; Cl-confidence interval. 
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Tab. 2 Genotype counts and p-values of the tested SNPs 


Genotype Counts 














Sample SNP ID Genotype P-value 
Cases Controls 
g.-847G>A GG/GA/AA 528/289/44 819/494/56 0.3046 
sf aA g.-888G>C GG/GC/CC T9A/141/5 1138/218/11 0.6012 
186951875 CC/CT/TT 766/169/5 1131/215/10 0.3495 
g.-847G>A GG/GA/AA 265/150/20 425/275/32 0.5707 
see g.-888G>C GG/GC/CC 409/70/2 614/113/5 0.7742 
186951875 CC/CT/TT 389/89/3 611/108/5 0.2340 
g.-847G>A GG/GA/AA 263/139/24 394/219/24 0.3341 
Prani ; g.-888G>C GG/GC/CC 385/71/3 524/105/6 0.7952 
186951875 CC/CT/TT 377/80/2 520/107/5 0.8258 





Genotypic p-values were calculated using PLINK under additive model (Purcell et al, 2007). 


significant results were observed (Tab. 1 and Tab. 2). 
Finally, we performed a meta-analysis by combining a 
previous study with a total sample size of 2 843 (Chen et 
al, 2002), but still failed to find a significant association 
between g.-888G>C and schizophrenia (OR=0.86, 95% 
CI=0.71-1.04, P=0.12) (Tab. 3). 

We further performed a haplotype-based analysis of 


the three SNPs with schizophrenia in all samples, female 


samples and male samples. No haplotypes were 
significantly associated with schizophrenia (Tab. 4). 
Notably, the global p-value of haplotype analyses in the 
female samples reached significance (P=0.007), which is 
likely due to a false positive effect caused by the low 


frequency haplotypes. 


Tab. 3 Meta-analysis of association between g.-888G>C and schizophrenia 


Numbers of C allele 


Numbers of G allele 








Study Ethnicity OR(95%CI) 
Cases Controls Cases Controls 
Chen et al, 2002 Chinese 38 50 520 460 0.67(0.43-1.04) 
Li et at, 2011 Chinese 151 240 1729 2494 0.91(0.73-1.12) 
Total 189 290 2249 2954 0.86(0.71—1.04) 


For the 2 association studies, fixed effect model based pooled OR=0.86, 


95% CI=0.71-1.04, Z=1.57, P=0.12; heterogeneity test: y’=1.45, df 











=1, P=0.23. 
Tab. 4 Analysis of the 3 SNPs for haplotypic association with schizophrenia 
Sample Haplotype fre_cases fre_controls CHISQ p-value Global p 
GCT 0.0811 0.0808 0.001 0.97 0.98 
All AGC 0.2215 0.2243 0.05 0.82 
GGC 0.6974 0.6949 0.03 0.86 
GCT 0.0771 0.0777 0.003 0.96 0.007 
Female AGC 0.2187 0.2321 0.59 0.44 
GGC 0.6840 0.6850 0.003 0.96 
GCT 0.0837 0.0840 0.001 0.97 0.93 
Male AGC 0.2199 0.2130 0.1485 0.70 
GGC 0.6965 0.7030 0.1063 0.74 


We performed a power calculation using the G*power 
program. The present sample size revealed >99% power to 
detect a significant association (a<0.05) with an effect 
size index of 0.1 (corresponding to a ‘weak’ gene effect). 


3 Discussion 


It has been reported that SNP g.-888G>C shows an 
association tendency with schizophrenia (P=0.08) (Chen 
et al, 2002). Our current study, however, showed no 
significant associations, suggesting SNPs in the promoter 
region of RELN confer no risk for schizophrenia or that 
there is a potential genetic divergence among regional 


Han Chinese populations. 

In our previous study, we observed multiple SNPs 
within RELN significantly associated with schizophrenia 
in Han Chinese. These SNPs were all located in the 
intron region and not causal SNPs (Li et al, 2011), 
however, suggesting the possibility of finding risk SNPs 
in the promoter region which had not been systematically 
screened before. We failed to find significant association 
between RELN promoter SNPs and schizophrenia, 
implying that the causal SNP was not located in the 
promoter region. However, there are several possibilities 
that could explain the negative results in this study. 
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Firstly, since we only sequenced 2 kb of the RELN 
promoter region, we could not detect other potential 
genetic variants risks not located in the -2 kb promoter 
region of RELN. Secondly, though our sample size was 
relatively large, it was smaller than the GWA study 
sample sizes and was unlikely to identify the risk variant 
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